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Abstract 

The PMNS neutrino mixing matrix C/pMNS is in general a product of two unitary ma- 
trices f/iep and Ui, arising from the diagonahzation of the charged lepton and neutrino 
^ , mass matrices, C/pMNS = U}^Uy. Assuming that Uy is a bimaximal mixing matrix, we 

Q I investigate the possible forms of U\ep. We identify three possible generic structures of C/iep, 

• which are compatible with the existing data on neutrino mixing. One corresponds to a 

hierarchical "CKM-like" matrix. In this case relatively large values of the solar neutrino 
mixing angle 6*801, and of |C/e3p = |(^>'^PMNs)e3pi are typically predicted, tan^ 0soi ~ 0.42, 
|t^e3p ~ 0.02, while the atmospheric neutrino mixing angle 0atm can deviate noticeably 
, from 7r/4, sin^ 20atm ~ 0.95. The second corresponds to one of the mixing angles in 

C/iep being equal to 7r/2, and predicts practically maximal atmospheric neutrino mixing 
sin^ 20atm — 1- Large atmospheric neutrino mixing, sin'^ 20atm ~ 0.95, is naturally pre- 
dicted by the third possible generic structure of f/iep, which corresponds to all three mixing 
angles in Ui^p being large. We focus especially on the case of CP-nonconservation, analyz- 
ing it in detail. We show how the CP-violating phases, arising from the diagonahzation 
of the neutrino and charged lepton mass matrices, contribute to the measured neutrino 
mixing observables. 
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1 Introduction 



After the spectacular results obtained in the experimental studies of neutrino oscillations 
in the last two and a half years or so (see, e.g., [TII2] for a summary), understanding the detailed 
structure and the origin of neutrino masses and mixing is of prime importance. Progress in 
understanding the origin of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix in 
the weak charged lepton current |HI can lead to a complete solution of the fundamental problem 
regarding the structure of the neutrino mass spectrum, which can be with normal or inverted 
hierarchy, or of quasi-degenerate type. It can also help gain significant insight on, or even 
answer the fundamental question of, CP-violation in the lepton sector. 

The fact that, according to the existing data, the atmospheric neutrino mixing angle is 
maximal, or close to maximal |3] , ^atm ~ vr/4, the solar neutrino mixing angle is relatively large 
12], ^0 ~ 7^/5.4 (sin^^Q = 0.30), and that the mixing angle 6 limited by the CHOOZ and 
Palo Verde experiments 0171 is small, sin^ 9 < 0.074 (99.73% C.L.) 2], suggests that the PMNS 
neutrino mixing matrix ^, f/pMNS, can originate from a 3x3 unitary mixing matrix having a 
bimaximal mixing form, f/bimax- In the case of CP-invariance in the lepton sector one has 
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The assumption of exact equality t/pMNS = f^bimax imphes 6q = 7r/4, which is ruled out at more 
than 5 s.d. by the existing solar neutrino data [3I2]- However, the deviation of ^0 from tt/4 can 
be described by a relatively small mixing parameter A = sin(7r/4 — it/ 5 A) = 0.20. Thus, 
the neutrino mixing matrix can have the form 

f^PMNS = Ul f/bimax ■ (2) 

It is natural to suppose that f/| and f/bimax in Eq- (0) arise from the diagonalization of the 
charged lepton and neutrino mass matrices, respectively. 

The existing data show also that |21ll] AttIq -C |Am^|, where Am^ > and Am\ are the 
neutrino mass squared differences driving the solar and atmospheric neutrino oscillations. This 
fact and the preceding considerations suggest further that f/bimax can arise from the diagonal- 
ization of a neutrino mass term of Majorana type having a specific symmetry. In the absence 
of charged lepton mixing {Ux = 1) this symmetry could correspond to the conservation of the 
non-standard lepton charge P L' = L^ — L^ — L^, where L/, I = e, fi, r, are the electron, muon 
and tauon lepton charges. The indicated symmetry cannot be exact: it has to be broken mildly 
by the neutrino mass term in order to ensure that Am^ ^ 0, and by the charged lepton mass 
term [Ux ^ 1), and/or by the neutrino mass term, in order to guarantee that ^0 7^ n/A. 

The possibility expressed by Eq. Q and described above has been discussed in the context 
of grand unified theories by many authors (see, e.g., jTUl El)- It has also been investigated 

^Throughout this article the case of three flavour neutrino mixing is considered. 
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phenomenologically first a long time ago (in a different context) in (H] and more recently in 
IT^ IT^ (see also J5j). Assuming that Eq. Q holds, we perform in the present article 
a systematic study of the possible forms of the matrix in Eq. which are compatible 
with the existing data on neutrino mixing and oscillations. The case of CP-nonconservation 
is of primary interest and is analyzed in detail. We show, in particular, how the CP-violating 
phases, arising from the diagonalization of the neutrino and charged lepton mass matrices, can 
influence the form of Ux and can contribute in the measured CP-conserving and CP-violating 
observables, related to neutrino mixing. The analysis presented here can be considered as a 
continuation of the earlier studies quoted, e.g., in HD 112 113 El ■ However, it overlaps little 
with them. 

2 The Deviations from Bimaximal Mixing 

We will employ in what follows the standard parametrization of the PMNS matrix: 



PMNS 



/ C12C13 S12C13 si3 ^^ 

-S12C23 - Ci2S23Sl3e"' C12C23 - 5125235136"^ 5230136"^ 
\ 512S23 - Ci2C23Si3e"^ -C12S23 - Si2C23Si3e"^ 0230136*'' / 



diag(l,e-,e^O, (3) 



where we have used the usual notations Cij = cos 6ij, Sij = sin^jj, 6 is the Dirac CP-violation 
phase, a and /3 are two possible Majorana CP-violation phases jlHllIIj. If we identify the 
two independent neutrino mass squared differences in this case, Am^i and Am'^i, with the 
neutrino mass squared differences which induce the solar and atmospheric neutrino oscillations, 
AtjIq = Am2i > 0, Am\ = Amg^^, one has: ^^12 = ^^soi, ^23 = ^atm, and 6*13 = 9. The ranges of 
values of the three neutrino mixing angles, which are allowed at 1 s.d. (3 s.d.) by the current 
solar and atmospheric neutrino data and by the data from the reactor antineutrino experiments 
CHOOZ and KamLAND, read 

0.35 (0.27) < tan2^,oi = tan^ ^12 < 0.52 (0.72) , 

|?7e3p = sin2^i3 < 0.029 (0.074) , (4) 
sin 2^atm = sin^ 2^23 > 0.95 (0.85) . 

As is well-known, the oscillations between flavour neutrinos depend on the Dirac phase S, 
but are insensitive to the Majorana CP-violating phases a and f3 [TBI H^!- Information about 
these phases can be obtained, in principle, in neutrinoless double beta decay ((/?/?) oi/-decay) 
experiments [201 1^11221 • 

Let us define the matrix U\ through Eq. ((21), f/pMNS = f/bimax, where the matrix t/bimax is 
given by Eq. (^Q). The latter would have coincided with the PMNS matrix if the neutrino mixing 
were exactly bimaximal. The bimaximal neutrino mixing can be obtained, e.g., by exploiting 
the flavor neutrino symmetry corresponding to the conservation of the non-standard lepton 
charge L' = (Lg — L^ — Lr) (see also [221 123)- The deviations of the neutrino mixing from 
the exact bimaximal mixing form can be described, as was shown in [8^, by a real parameter 
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A ~ 0.2, which was introduced in the following "flexible" parametrization of three elements of 
f^PMNs: 

f/e2 = /|(1-A) , ?7e3 = AA", U = {I - B X^^) (5) 

where A and B are real parameters of order one. The integer numbers m and n can be chosen 
according to the improved limits on, or precise values of, sin^ 6'i3 and sin^2^^23- The fact that 
solar neutrino mixing is now conflrmed (at more than 5 s.d.) to be non-maximal means that 
A 7^ and the best-flt value of tan^ 6*801 = 0.43 corresponds to A ~ 0.23. 

Taking the purely phenomenological parametrization ^ at face value and expressing the 
PMNS matrix in terms of A, A and B, one can solve Eq. (0) for U\ in order to obtain the 
physical ranges of values of the parameters involved. It was found in P] that the structure of 
U\ can at leading order be the unit matrix plus corrections of order A. 

In the present article we perform a systematic study of the possible forms of U\ allowed 
by the current data. The neutrino mixing phenomenology can be described, for example, by 
U\ having a "CKM-like" structure, as has been discussed in ^3 ^3 |H]. Here we show, in 
particular, that the restriction to a hierarchical form for U\ is not necessary. Counter-intuitive 
forms of [/a, and/or the presence of CP-violating phases, can produce naturally the observed 
deviation from 7r/4 of ^soi as well as the requisite smallness of | Uez \ ■ 

3 The Case of CP Conservation 

In this case U\ is an orthogonal matrix. We will use the standard parametrization (see 
Eq. (jni) for U\. Let us denote the angles in U\ as 6'i2i 6*23, and 6[^, and deflne sin^-^- = Xij. 
We shall treat A12, A23 and A13 as free parameters without assuming any hierarchy relation 
between them. In the case of CP-conservation under discussion one can limit the analysis 
to the case < 6[j < vr, and correspondingly to < \ij < 1. The results for Xij < can 
be obtained formally from those derived for Aj^ > by making the change i) A23 —>■ — A23, 
and/or ii) A13 —X13 as long as one keeps A12 > 0. The change A12 — A12 should be 
done simultaneously with the change A13 — A13 and in this case only the "solutions" with 
1^13 1 > IA12I should be considered. In the latter case both signs of A23 are possible. 

Using Eqs. ^ and (jH), we have determined the regions of values of A12, A23 and A13 for 
which the values of the neutrino mixing angles 6*12 = 6*801, ^23 = ^atm, and 6*13 = 9, lying within 
their la and 3cr ranges, Eq. (^, can be reproduced. The results are shown graphically in Fig. 
n One can clearly identify three different cases: 

(i) all Xij < 0.35 ("small"), 

(ii) A23 = +1 and Ai2,i3 ^ 0.35 ("small"), 

(iii) all Xij > 0.40 ("large"). 

We shall discuss next these three cases in detail. The resulting structures of Ux will turn out 
to be rather different in the three cases. 
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3.1 Small X,j 

Multiplying Ul by f/bimax and expressing tan^ 6'soi, l^^esP and sin^ 2^atm in terms of the three 
small parameters Xij, one finds: 



tan^ 9. 



sol 



\u, 



e2\ 



l-2v^(Ai2 

A12 + A 



-^13) + 4 (A12 — Ai3 



U, 



eSl 



sin 26'a,tTn = 



13 _ A^a/^ 

4|[/„,,P|f/,3|' 



2-2v^A23(Ai2 + Ai3)+0(A3) 

12 - A13) + 0(A3) 



V3| 

:i - If/, 



e3| 



(6) 

The key quantities are therefore A23 and the sum and the difference of A 12 and A13. To leading 
order in the A-parameters, the deviation from maximal solar neutrino mixing and zero t/es 
are proportional to A, while the atmospheric neutrino mixing is close to maximal with only 
quadratic corrections to (1 — sin^2^atm)- The requirement of sin^ 26'atm > 0.85 (0.95) leads to 
order A^ to the restriction A23 ~ 0.19 (0.11). 

Assuming that the oscillation parameters lie within their Icr allowed ranges and neglecting 
terms of order A^, one finds from Eq. ® that the term (A12 — A13) has to be positive and 
rather large (~ 0.2) in order to ensure the deviation of ^^soi from it /A. At the same time 
(A12 + A13) has to be smaller than ~ 0.2 in order to satisfy the limit on \Ue3\- These two 
conditions imply (for Ai2,i3 > 0) a hierarchy of the form A12 3> A13. The remaining parameter 
is A23 ^ (A12 + Ai3)/2 ~ A12/2. 

In Fig. 121 we show scatter plots of the values of Xij, obtained by requiring that the corre- 
sponding values of neutrino mixing parameters tan^ 6'soi, sin^ 2^atm and |?7e3p lie within their 
la allowed ranges, Eq. Q. One can see from the left panels in Fig. |21that A12 — (0.21 — 0.26), 
A23 ~ 0.10 and A13 ^ 0.03. Thus, one has A23 ^ A12/2 and A13 ^ A12/2. In Fig. |2lwe also 
display the resulting correlations between the neutrino mixing parameters for which the la 
allowed ranges were used. From the panels in the right column one sees that tan^ ^soi ~ 0.42 
and |f/e3p ~ 0.017. Most of the points tend to lie at rather large values of tan^ 6'soi and |f/e3p- 
Taking the 3a range in Eq. leads to the lower limits of tan^ 6'soi ~ 0.35 and of |f/e3p ~ 0.003. 

We consider next two representative (and to a certain degree typical) forms of U\. A viable 
possibility is the existence of "hierarchical" relations between A12, A23 and A13. One can have, 
for instance, A23 ~ XI2 -^is ~ -^125 "with A12 = A and, e.g., A ~ (0.20 — 0.30). This implies 
the following form of U\: 



( 1 - AV2 



V A3 



A A3 \ 
I-AV2 A2 
-A2 1 / 



C(A^ 



(7) 



The structure of Vx is close to that of a diagonal matrix, and is similar to the structure of the 
CKM matrix. Given the hierarchy in the charged lepton masses, the CKM-like form of [/a, Eq. 
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(|7j), is rather natural. One has in this case: 

tan^^soi ^ 1 - 2v^A + - 2v^A3 



\U, 



e3 



72' 

sin2 2^,tm^ 1-4A4 

plus terms of higher order in A. Consequently, this scenario "predicts" solar neutrino mixing 
and |f/e3p close to their currently allowed maximal values, and atmospheric neutrino mixing 
close to maximal. 

More specifically, the hierarchy of the A-parameters we are considering, A12 ^ A23,i3, leads, 
as it follows from Eqs. (jHl) and (jH)), to the interesting correlation: 

tan^ 9,oi -l-4\Ue3\ + S\Ue3\^ -SlUesl"" . (9) 

The upper limit of tan^ 6'soi < 0.52 (0.72) implies in this case a lower limit on |t/e3p ~ 
0.027 (0.006) ^. This explains qualitatively the lower limit on |f/e3p as observed in Fig. |21 
and mentioned earlier. Another correlation related to the CKM-like form of Ux, Eq. ((7j), reads: 

sin2 20,tm^l-16|f/e3|'. (10) 
Thus, the deviations from maximal atmospheric neutrino mixing are determined in the case of 
the hierarchical relations A12 = \ ^ A23,i3, A23 ~ A^, by the magnitude of |?7e3|^- For, e.g., 
|[/e3p = 0.029 (0.05), we have sin2 2^atm ^ 0.975 (0.96). 

The alternative possibility is that of a mild "hierarchy" between A12 and A23. We can have, 
for instance, A23 — A12/2, A12 = A. Taking also A13 ~ A^ ^ A12, one finds: 

tan^^soi^ l-2V2A + (4- V2)A2-2V2A=^ , 

\U.\ ^ ^ . (11) 
sin2 2^,tm^l-A2~l-2|[/e3p . 

The atmospheric neutrino mixing angle ^atm can deviate more from 7r/4 than in the hierarchical 
case: for |?7e3p = 0.05 we have now sin^ 26'atm — 0.90. This provides a possibility to distinguish 
the case of A23 — A12/2, Ai2,23 ^ A13, from the hierarchical one A13 ~ Afg < A23 ~ A^2 < -^12 con- 
sidered above. For A ~ 0.24 one finds from Eq. (jllj) tan^ ^soi — 0.43, which is quite encouraging 
since values of A ~ 0.24 could be interpreted in terms of ^Jm^jm^ ~ a/O. 105/1. 77 ~ 0.24. 

3.2 The Case of A23 = 1 

If A23 = 1 and Ai2,i3 are "small", Ai2,i3 ~ 0.35, we find: 

tan2 1 - 2v^ (A12 + A13) + 4 (A12 + \xzf + ^(A^) , 

f/,3 ^ ^^^^^ + 0[}?) , (12) 
sin2 2^atm ^ 1 - i (A?2 + A^g)^ + A12 A13 (A?2 - A^g) + ^(A^) . 



^That the upper limit on tan^ 6*801 implies a significant lower limit on |C/e3p in the case of CKM-like form of 
\J\ was noticed also in 
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The correlations of A12 and A13 as well as of the mixing parameters for A23 = 1 are shown 
in Fig. El We used the la allowed ranges of the oscillation parameters, Eq. (jlj). As Fig. El 
shows, sin^ 26'atm is practically 1 in this case, sin^ 2^atm ^ 0.9965. The lower limit is reached 
for relatively large |t/e3p — 0.025 and relatively small tan^ ^soi — 0.32. However, such small 
deviations of sin^ 2^atm from 1 are extremely difficult to measure ^. For A23 = 1 and "small" 
-^12,13, U\ has the following structure: 
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v 



\ (-^12 + -^13) 
— Ai3 

A12 



A12 

— A12 Ai3 
-1 + U?2 
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13 



2 ^13 




+ 0(A' 



(13) 



All three possibilities A12 S> A13, A13 ^ A12 and A12 ~ A13 are allowed. In any of these three 
cases the 23 submatrix of \J\ displays an antidiagonal structure, as seen in Eq. ()13p . 



3.3 Large A^j 

It is also possible that all Ajj = sin^^^ are relatively large, Ajj ^ 0.40. From Fig. HI one 
sees that in this case A12 ^ (0.92 - 0.99) [0.85 - 0.99], A13 ~ (0.57 - 0.77) [0.47 - 0.85], and 
A23 — (0.57 — 0.92) [0.40 — 0.99], when the 1 s.d. [3 s.d.] ranges of the oscillation parameters 
are used. It turns out that, cos 6*23 > 0, whereas cos6'']^2,i3 can take both signs. 

We limit our discussion to the case of all A,, being "large": \ij > 1/^2. We can use the 
deviations of A12 from 1, and the deviations of A1323 from as small parameters to get 

convenient expressions for tan^6'soi, |?7e3p and sin^26'atm: 

A12 = 1 ~ ^12 

Ai3 = 73 (1 + 'rs) (14) 
A23 = 75 (1 + ^23) ' 

where ef^ ^ 0.08 (0.15), ej^ < 0.30 (0.40) and e?3 < 0.09 (0.20), where the 1 s.d. (3 s.d.) ranges 
of the oscillation parameters were used. One obtains for sufficiently small eij: 

tan2 ~ 1 _ 4 ei2 + 8 + 0{e^) , 

Ue3 ^ 4, - ei2 + Oie^) , (15) 

sin2 2eatm=^l + 0(e^) . 

The terms of higher order in e we do not give can be sizable. The scatter plots do not show 
any new important correlation. In particular, sin^ 2^atm — 0.95 is possible also in the case 
under discussion. Both €23 and €13 can be zero, whereas ei2 cannot. No hierarchy is implied. If 

ei2 > ei3,23, we get {Uesl - ^12(1 - e?2/4) + C'(e^), and 

tan^ 9,ol ~ 1 - 4 |f/e3| + 8 If/esp - 16 |f/e3|' . (16) 
■^The requisite precision may be achieved in experiments at neutrino factories |25| . 
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For all Xij "large" , Xij > 1/ V^, U\ has the form: 



ei2 



1 

V2 



f2 -|_f2 



-1 £12 



''12^^ '=23 

"=23 '=12 — 1 
2 



"T + ^12 + I (e?2 - ^13 ~ e 

2 2 

_i tl2. I '=23~'=12 _ 

\ V2 V2 2 '=12 -T 2 l'=12 '^IS 

All entries in U\, except the 11 entry, are large. 



23) 



23 y 



+ 



^(£23 -e 



^=23 



13; 



+ 0(e3) . (17) 



+ 4) / 



Let us summarize the results obtained under the assumption of CP-conservation: 

• If all Xij are small, one typically has tan^ ^soi ~ 0.42 and |t/e3p ~ 0.02. The matrix U\ 
has a "CKM-like" structure. One can have sin^ 2^atm — 0.95, i.e., sin^ 26'atm can deviate 
noticeably from 1. 

• If A23 = 1 and Ai2,i3 are small, one practically has sin^ 26'atm — 1, the possible devia- 
tions from 1 being exceedingly small. Any deviation of sin^ 26'atm from 1 would imply 
rather large |f/e3p and relatively small tan^6'soi- The 23 block of the matrix U\ has an 
antidiagonal (quasi-Dirac like) form. 

• If all Xij are large, Xij ^ l/"\/2, the matrix U\ displays the unusual structure given in Eq. 
()17|). with all elements, except the 11 element, being large. Also in this case we can have 
sin^ 2e^tm ^ 0.95. 



3.4 Implied Forms of the Charged Lepton Mass Matrix 

The matrix Ux arises from the diagonalization of the charged lepton mass matrix miep. 
Having the approximate form of Ux and using the knowledge of the charged lepton masses me, 
and m^, we can derive the structure of two different matrices related to miep. The first one 
is the matrix miep "^lep? which is diagonalized by Ux, i.e.. 



"^lep "^Lp = Ux (mJpS)2 Ux^ 



where tti^I^ is a diagonal mass matrix containing the charged lepton masses me, m^ and m,-. 
The second matrix is 

(19) 



mg^ Ux 



T 



Under the assumption of miep being symmetric, which is realized in some GUT theories, the 
matrix in Eq. (fTIHl coincides with the charged lepton mass matrix. 



Consider first the matrix miep ^ 



t 

lep' 



In the case of "small" Xij, we find that 



/ ml + ml A?3 + ml Xj^ ml A13 A23 + ml Xu ml X 



ml + ml X^ 



23 



V13 



m^l, A12 Ai3 \ 



m^A 



23 



m 



^ Ai5 Ai 

^ -^12 '^IS 



(20) 



m: 
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plus terms of order A^. Thus, miep mj,p is close to a diagonal matrix. In the special case of 



A 



13 



A23 ~ we get: 



/ ml + w? \^ mi X \ 



mf, A 



mi 
V / 



(21) 



The zero entries can be terms of order m^. The matrix miep m\^^ has clearly a hierarchical 
structure. 



When A23 = 1 and Ai2,i3 are "small", one obtains: 



/ ml + ml A?3 + ml A?2 ml A13 



-m^ A12 \ 



A 



13 



-ml A12 



A12 Ai3 



A12 Ai3 



(22) 



In contrast to the previous case, now the 22 element is equal to ml and the 33 element is given 
by m^. The 12 block of miep mjep is hierarchical, while the 23 block has "inverted hierarchy-like" 
form. 

Finally, if all Xij are "large", we can expect that micp m\^^ is to a good approximation a 
"democratic" matrix. Indeed, one finds: 



t ^7 



I 1 1/v^ 1/^2 \ 

I/V2 1/2 1/2 
\ I/V2 1/2 1/2 / 



(23) 



plus terms of order e. 



From the definition of the two characteristic mass matrices in Eqs. fjl8|) and (fT^ it follows 
that the expressions for Ux m^l^ Uy^ in the three cases we are discussing can be obtained from 
the corresponding expressions of miep "^lep given above by replacing mf with m/ in the latter. 



4 CP Nonconservation 



4.1 General Considerations 

It is well known that if the neutrinos with definite mass are Majorana particles, the PMNS 
matrix contains 6 physical parameters, three mixing angles, and three CP-violating phases. 
The latter can be divided in one "Dirac phase", 5, and two "Majorana phases", a and (3. The 
Majorana phases appear only in amplitudes describing lepton number violating processes in 
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which the total lepton charge changes by two units. In general, as we have already discussed, 
one has 



f/pMNS = t/ilpf/. • (24) 

where Uiep and U^, two 3x3 unitary matrices: f/iep arises from the diagonalization of the charged 
lepton mass matrix, while Uu diagonalizes the neutrino Majorana mass term. Any 3x3 unitary 
matrix contains 3 moduli and 6 phases and can be written as |2nj 

U = e'^PUQ, (25) 

where P = diag(l, e"^^, e*'^^) and Q = diag(l, e*^^, e*''^) are diagonal phase matrices having 2 
phases each, and f/ is a unitary "CKM-like" matrix containing 1 phase and 3 angles. The 
charged lepton Dirac mass term, miep, is diagonalized by a bi-unitary transformation: 

miep = ^Lmg^^i, (26) 

where Ul,r are 3x3 unitary matrices and mf^^^ is the diagonal matrix containing the masses of 
the charged leptons. Casting Ul,r in the form ((23), i.e., Ul,r = e**-^'^ Pl,rUl,rQl,r, we find 

miep = e^(*-*-) Q, U, P, mf^ Q^^ Ul pt . (27) 

The term P^ '^fep^ contains only 2 relative phases, which can be associated with the right- 
handed charged lepton fields. The three independent phases in e*^*^~*^^(5L can be absorbed 
by a redefinition of the left-handed charged lepton fields. Therefore, f/iep is effectively given by 
Ul and contains three angles and one phase. 

The neutrino mass matrix m^, is diagonalized via 



U,,m'!!''^UT . (28) 



The unitary matrix can be written in the form ()25|1 . It is not possible to absorb phases in 
the neutrino fields since the neutrino mass term is of Majorana type fHlEI- Thus, 

t/pMNs = f//ep = e'"" f/ip P. f/. . (29) 

The common phase $jy has no physical meaning and we will ignore it. Consequently, in the 
most general case, the elements of ?7pmns given by Eq. are expressed in terms of six real 
parameters and six phases in [/lep and Ui, ^ . Only six combinations of those — the three angles 
and the three phases of f/pMNS; are observable, in principle, at low energies. Note that the two 
phases in Qy are "Majorana-like" , since they will not appear in the probabilities describing the 
flavour neutrino oscillations [^E]. Note also that if f/iep = 1, the phases in the matrix Py 
can be eliminated by a redefinition of the charged lepton fields. 

^Using a different representation of the unitary matrices C/jep and J7^, the authors of ^1] came to the 
conclusion that t/pMNS is expressed in terms of six real parameters and seven phases of C/iep and Uy. Writing 
the unitary matrices in the form of Eq. (|25|l allows to reduce the number of phases by one. 
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The requirement of bimaximality of If^, implies that Ui, is real and given by Eq. (^J). In 
this case the three angles and the Dirac phase in the PMNS matrix t/pMNS will depend in a 
complicated manner on the three angles and the phase in U\ep and on the two phases in P^. 
The two Majorana phases will depend in addition on the parameters in Q^,. 

It should be emphasized that the form of ?7pmns given in Eq. is the most general one. A 
specific model in the framework of which the bimaximal 11^ is obtained, might imply symmetries 
or textures, along the lines of, e.g., [32], in m^, which will reduce the number of independent 
parameters in U^. We will comment in Section on such possibilities. 

In the scheme with three massive Majorana neutrinos under discussion there exist three 
rephasing invariants related to the three CP-violating phases in Upmns, ^, « and /3 |2Zll2Ell2ni 
ISniE]- The first is the standard Dirac one Jcp [27J, associated with the Dirac phase S: 

JcP = lm{UeiU^2U:^U;,} . (30) 

It determines the magnitude of CP-violation effects in neutrino oscillations j2Hl- Let us note 
that if Uiep = 1 and Ui, is a real matrix, one has Jcp = 0. 

The two additional invariants. Si and 5*2, whose existence is related to the Majorana nature 
of massive neutrinos, i.e., to the phases a and P, can be chosen as [211 EI] (see also [21]) ^: 

^1 = Im {f/,1 U:,} , ^2 = Im {f/e2 U:,} . (31) 

If Si 7^ and/or 5*2 7^ 0, CP is not conserved due to the Majorana phases P and/or (a — /?). 
The effective Majorana mass in (/?/3)oi/-decay, |<m>|, depends, in general, on Si and S2 PT] 
and not on Jcp. Let us note, however, even if S'1^2 = (which can take place if, e.g., |?7e3| = 0), 
the two Majorana phases a and (3 can still be a source of CP-nonconservation in the lepton 
sector provided Im {Uei Ul^} ^ and Im { f/^2 f/^s } 7^ [SH IS2] • 

Let us denote the phase in [/lep by if). We will include it in ?7iep in the same way this is done 
for the phase 5 in Eq. We will also use Py = diag(l, e*''^, e**^) and = diag(l, e*'', e*'^). 
This means that the Dirac phase 6, which has observable consequences in neutrino oscillation 
experiments, is determined only by the phases ip, and u. The Majorana phases in ?7pmns, « 
and /3, receive contributions also from the two remaining phases p and a. Allowing the phases 
(5, a and j3 to vary between and 2-7?, permits to constrain (without loss of generality) the 
mixing angles 9ij to lie between and tx/2. In Fig. Hjwe show the result of a random search for 
the allowed regions of the values of the parameters The regions present in Fig. ^ are again 
allowed. Clustering of points around these regions is also observed. The effect of the phases is 
that values of the parameters, located in areas between the regions corresponding to the case of 
CP-conservation are allowed. Let us consider the physical observables of interest in the three 
main regions of the parameter space we have identified earlier. 

^We assume that the fields of massive Majorana neutrinos satisfy Majorana conditions which do not contain 
phase factors. 
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4.2 Small X,j 

Let us choose A12 = A, A23 = AX and A13 = B \ with A, B real and of order one ^. In this 
case tan^ 6*801, t^es and sin^ 26'atm are given by 

tan^^soi^ l-2v^(c^-5c^_^)A + 2(2(c^-5c^_^)2- V2A(c^ + 5c^_^,)) A' , 

. (a - .l'^^^^:-^ A^) . (32) 

sin^ 2^atm ^ 1 - 4 ^2 c^_^ A2 

plus terms of order O(A^). We introduced the obvious notation = coscf), s^_^ = sin(a; — ip), 
etc. Note that the possibility of all three Xij being of the same order was not possible in case of 
conservation of CP. The rephasing invariant Jcp, which controls the magnitude of CP-violating 
effects in neutrino oscillations, has the form: 

Jcp ^ ^ is<t, + B + A (s^_20 + B S2^-^-^) A) + 0{X^) . (33) 

Thus, not only the phase of Uiep, but also the phases uo and from contribute to Jcp- 
Using Eqs. (fH^ and (jHSI) it is not difficult to convince oneself that the magnitude of Jcp is 
controlled by the magnitude of \Uez\- Indeed, if, e.g., c^;-^-^ = — 1 and S = 1, so that the term 
~ A in the expression for Uez vanishes and to leading order [/e3 ~ —\/2 Ac^_^ A^, the term ~ A 
in Jcp also vanishes and we have Jcp ^ Ues- 

The invariants Si and 5*2 which (for |?7e3| 7^ 0) determine the magnitude of the effects of 
CP-violation associated with the Majorana nature of massive neutrinos, read: 

Si ^ h (s^+cr + B 

A ^^^^ 
S2 — {s4>~p+cT + B Stj_^_p+o-) , 

where we gave only the terms of order A. We see that all five phases in JJiep, Pu and Q^, and 
not only the phases p and a in Q^^, contribute to Si and S2- In the leading order expressions 
for 5*1^2, Eq- (jS3), the five phases enter in three independent combinations. Moreover, in the 
case of p = 0" = 0, which is realized if the bimaximal mixing structure of is associated (in 
the limit of U\ep = 1) with the approximate conservation of L' = Le — — (see Subsection 
14. 6p . we have to leading order in A: 

Thus, the magnitude of the CP-violating effects in neutrino oscillations is directly related in 
this case to the magnitude of the CP-violating effects associated with the Majorana nature of 
neutrinos. To leading order, both types of CP-violating effects are due to two phases, and 
{uo — ip). We will consider a model in which p = o" = in Subsection 14.61 

^Hierarchical Ay correspond in this paramctrization to sufficiently small A and/or B. 
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In Fig. we show the correlations between the observables tan^6'soi, |f^e3p and sin^2^atm5 
when a hierarchy between the three A parameters of the form we have considered in Subsection 
13 .21 Ai2 = A, Ai3 = A'^ and A23 = A^, is assumed. The la allowed ranges of tan^6'soi, l^esP and 
sin^ 26'atm5 were used to constrain A and the CP-violating phases. As it follows from Fig. 
tan^^soi can take values in the whole la interval allowed by the data, Eq. (0]). This is in sharp 
contrast to the case of CP-conservation, in which the correlation between the values {Ue^l and 
the deviation of tan^ 6^801 from 1 leads, as Fig. El shows, to the lower limit tan^ ^soi ~ 0.48. We 
find also that if CP is not conserved, the following lower bound holds: |t/e3p ~ 0.018. There 
are interesting correlations between tan^ 6*801, If^esP and sin^ 26'atm- For instance, tan^ 9so\ < 0.40 
implies |[/e3p Z 0.026 and sin^ 2^atm ^ 0.975. 



4.3 The Case of A23 = 1 

This case corresponds to "small" Ai2,i3, Ai2,i3 ^ 0.35. Introducing A12 = A and A13 = B\, 
with B real, we find: 

tan^ e,oi -l-2V2{c^ + B c<^_^) X + A{c^ + B c^_^f A^ , 



?7e3^ A-y^(l + 52-25c^„^+^) + 0(A3) , (36) 

sin^ 2^atm -1-\{1- B^ -2B c^^^^^f A^ . 

The formulae for the mixing parameters can actually be obtained from the corresponding 
formulae for all Ajj small, Eq. (j32|) . by setting A = and making the change B ^ —B, (p ^ u. 
The rephasing invariant Jcp reads: 

Jcp -^{s^-B s^.^) + 0{X^) . (37) 

The Majorana counterparts to Jcp-, Si and 5*2, are given by 

' (38) 



Si — — ^ (st^+CT — B S0„^+ct) 



S2 — 2 {^^^~p+cr B S(j,^^-p-^cr) ■ 

For p = a = one finds that to order A the relation Si = S2 = —2^/2 Jcp is valid. Relatively 
small uj is favored in this case. Furthermore, atmospheric neutrino mixing is predicted to be 
close to maximal, to be more precise, the lower limit sin^ 26'atm ~ 0.97 holds. 



4.4 Large \ij 

We consider Xij > 1/ \/2 and choose again as small expansion parameters the e^j introduced 
in Eq. (|T^ . For the oscillation observables one finds 

tan^ ^3oi ^1-4 + ^12 + 2 ^ + ^ ''^^ ^2 , 

cos(cu - 0)/2 cos\uj - 0)/2 ^2 ' 

f/e3 ~ sin(c<j — 0)/2 + ei2 cos(c<j — 0)/2 , (39) 
sin2 2^atm^l + 0(e') . 
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The atmospheric neutrino mixing is again very close to maximal. The term ~ sin(co' — </))/2 
in Ue3 is not suppressed by positive powers of eij. This suggests that (a; — 0)/2 should be 
relatively small. For u = (p the expressions for the three rephasing invariants, associated with 
CP-nonconservation, take a rather simple form: 

^cP--|ei, + ^6^3, (40) 



and 

V2 V2 



In the case of a = p = we have 5*1 ~ — 5*2, but there is no simple relation between 5'i_2 and 
Jcp- 

Note that in all three cases under study the connection between Si, S2 and Jcp is different. 

4.5 Special Cases 

It is instructive to examine the interplay of parameters when some of the three phases in 
f/iep and Pj, are zero. It turns out that there are no drastic consequences in the cases of "small" 
hierarchical Xij, of A23 = 1, and when all Xij are "large". Interesting correlations happen, 
however, in the case of "small" when the three Xij are of the same order. To illustrate 
this we choose again A12 = A, A23 = AX and A13 = B X with A, B real and of order one. We 
fix A = ^Jm^^jm^ ~ 0.24 and vary A,B between 0.4 and 1/-\/A. In doing so we observe that 
values of A and B slightly below 0.5 are favored, indicating a mild hierarchy in A^. Inspecting 
Eq. (j32j) one finds that the dependence on lo of the oscillation parameters is rather weak. If 
uj = one finds from Eq. that tan^ 6*801 will be too large, tan^ 6*801 ~ 0.58. We show in 
Fig. inithe scatter plots of the correlations between tan^^soi and |t/e3p, and between |?7e3p and 
JcPi which turn out to be the most interesting ones. The \a allowed ranges of values of the 
neutrino mixing parameters given in Eq. (0)) were used. If the charged lepton mass Lagrangian 
conserves CP, i.e., if = 0, there exists a lower bound on |f/e3p ~ 0.005. In this case |t/cp| can 
be relatively large, for instance, \ Jcp\ — 0.035 for |f/e3p = 0.025, but is also allowed to be zero. 

If CP is conserved by the neutrino mass term, i.e., if = = 0, a simple correlation 
between tan^ 6*301 and |f/e3p exists again: smaller |t/e3p implies smaller tan^ 6*301 (see Fig. IHl). A 
similar limit on |t/e3p as in the case of = applies. Most remarkably, the area of allowed 
values of Jcp "bifurcates" as a function of |f/e3p and for values of |?7e3p ~ 0.008, Jcp = is no 
longer allowed. For |f/e3p > 0.010, for instance, one has \ Jcp\ ~ 0.01. This can be understood 
as follows. If uj = (f) = holds, the requirement of vanishing Jcp (see Eq. fj33|) ) is fulfilled for 
ip = 0. Then, one sees from Eq. (fH^ that to order A, the equality 

tan^^soi^ 1-4^— ||f/e3| (42) 
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holds. For the values of B considered, B = (0.40 — 2.0), this equality is incompatible with the 
la upper bounds on tan^ 6*801 and |^7e3p- 

The "bifurcation" happens also when (p = ip = holds (Fig. IH)). The lower limit on |?7e3p 
in this case is roughly by a factor of 4 smaller than in the preceding ones. Though not shown 
here, there can also be interesting constraints on sin^ 2^atm when some of the phases vanish. If, 
e.g., u = (p = holds, one sees from Eq. (j22I) that 1 — sin^ 26'atm receives a correction of order 

~ 0.06. It can be shown that typically sin^ 26'atm ~ 0.98 in this case. 



4.6 Specific Models 

We shall consider next specific models in which Eq. (I24j) arises. We will also investigate 
how the CP- violating phases appear in the matrix Uu- 

It is well-known that when ?7iep = 1, bimaximal neutrino mixing can be generated by the 
following neutrino mass matrix ES] 



m 

7! 



/o 


1 


1 \ 


1 








V 1 





0/ 



(43) 



The latter has eigenvalues of and ±m. This matrix has a flavour symmetry which corresponds 
to the conservations of the lepton charge V = Le — — Lr. One can lift the degeneracy of 
the two mass eigenvalues by adding a small perturbation em/\/2 in the /xr entry. This yields 
a mass matrix known from the Zee model j33j (see also, e.g., (SHI EHl EZI EH])- 

It is well-known that there is no CP-violation in the Zee model since one can absorb all 
three possible phases in the charged lepton flelds. As we will see, this is true only when the 
charged lepton mass matrix is diagonal. One finds that in this case Am\ ~ (1 — e/ \/2) and 
Amg ~ \/2rn? e. This fixes 



Am^ 1 7-10-^eV2 



v/2 Ami V2 2 ■ lO'^ eV^ 



0.025 , (44) 



where we used the best-fit values found in the recent analyzes in [21 E]. For e 7^ the solar 
mixing is slightly reduced from its maximal value jSHlEZl: we have tan^ ^soi — 1 — e/ \/2 ^ 0.98. 
This is incompatible even with the 5(j allowed range of values of tan^ 6^301 (see, e.g., |2j). 

In order to examine the implications of existence of phases in my, we will assume that the 
three non-vanishing entries in m^ are complex and will denote the phases of these entries — 
the e/i, er and /ir, by a', [3' and 7', respectively. The neutrino mass Lagrangian then has the 
form: 

£ = (iT^z/^^e-"' +77^i.;:^e-^' + eI7;Z^^;:^e-^') + h.c. (45) 



^Note that, as is shown in Fig.|Hl one has Jcp = for the minimal aUowed value of |C/e3p — 0.006. As can 
be seen from Eq. H33|) . the latter is reached in the case under discussion for the CP-conserving value of ■0 = tt, 
for which Jcp vanishes. 
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where ufff = C Vei^ ■, C being the charge conjugation matrix. We can redefine the fiavor neutrino 
fields as follows: = e*"'i/;^, / = e,fi,T. Choosing 

ae = -i(«' + /?'-7') , 

«M = -|(«'-/^' + 7') , (46) 
ar = -i(/3'-a' + y) , 

we get in terms of the fields z/^^: 

^=-^{K~Ar + K~l <r + e^L<R)+i^-0- (47) 

Now the neutrino mass matrix contains only real entries and is diagonalized by an orthogonal 
transformation with an orthogonal real matrix O. The matrix O has to a good approximation 
the bimaximal mixing form, Eq. (^). The mass eigenstate Majorana fields Xj> j = 1)2,3, are 
related to the fields u'^ = (i/^^, z/^^, u'^^^f via v'^ = Oxl, where xl = (Xil, X2L, XslV- Since u'^ 
can be written as Pj ul, where ul = (j^eL, J^^jL, ^tl)^ contains the original fiavor neutrino fields, 
and Pi, = diag(e*"% e*"'', e*"^), we have = P^O xl- We can factor out one of the three phases 
in Py, i.e., 

= e-t("'+/^'-V) diag(l, e-^(^'-^'), e'^^^'""')) = e"*"^ P, . (48) 

Thus, comparing the resulting formula with Eq. (j29j) . we see that in this special case of only 
three phases in m^, for the phase matrix holds: = 1. Consequently, the characteristic 
leading order relations between Si, S2 and Jcp as discussed in Sections 14.21 to lOl will hold in 
the three generic structures of Ux discussed above. In terms of the notations used earlier in 
this Section we have = (/?' — 7') and u = {a' — 7')- Thus, there are altogether 3 CP-violating 
phases in f/iep and in the specific model we are considering. In the case of a' = P', one has 

(f) = UJ. 

The mixing parameters tan^ 6'soi, Ues and sin^ 26'atm are functions of the 3 angles and 1 phase 
in U\ep and of the 2 phases in P^. The formulae and results we have derived in Subsections 
14.21 to 14.51 are valid in the specific model under discussion if one sets p = a = and takes 
into account the fact that (p = {(3' — 7') and uj = {a' — 7'). If f/icp = 1, the phases in P^ are 
unphysical and can be absorbed by the charged lepton fields, CP is conserved in the lepton 
sector and Jcp = Si = S2 = 

The neutrino mass term, Eq. ()45|) . leads to a neutrino mass spectrum with inverted hierar- 
chy: m2 — —mi, -C m2, and Am^ = m| — ml > 0, Am\ = m^ — m| ~ m^ S> Ami,. The 
effective Majorana mass measured in (/3/3)oi/-decay, |<m>|, is given approximately by (see, 
e.g., EH) 

V2m\ 1 1 - P e*^"'-'^'-^) I , case I 

(49) 

v^mA |l + Pe*(^'-"'-'^)| , case II 

where case I and case II refer respectively to A12 = A, A23 = A\, A13 = P A, and to A12 = A, 
A23 = 1, Ai3 = P A. Note that |<m>| does not vanish in the limit of e = 9J. Comparing the 



\<m>\ 



Am\ 



U^i - u. 



e2 
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expressions for |<m>| with that for JcPi 

{ Tm^^P'-y + ^ Sq/^y^^) a , case I , 
Jcp ~ <^ (50) 
I 471 (-^"'-7' ~ ^ sp>_Y_^) A , case II , 

one finds that all three phases (a'— 7'), — 7') (or (p, uj) and ip contribute to both, the "Dirac" 
phase 6 and to the single Majorana phase that enters into the expression for |< m>| in the case 
of neutrino mass spectrum with inverted hierarchy fl^ \'2'2\. 

If we set f3' = a', only ip contributes to |<m>|, while Jcp is a non-trivial function of the two 
phases (a' — 7') = (p and {a' — •y' — ip) = {(p — tp). If, however, /?' = 7', then = 0, = (a' — /?'), 
and we have to leading order in A in case I: Jcp — \ B Sq,'-/3'-^/(4-\/2) ~ Si/ (2-\/2) ~ S2/ (2\/2). 
The effective Majorana mass as given in Eq. (|49|) can be written as: 

|<m>| ~ m |cos26'soi — ^ 8 J^pI . (51) 

The same relation holds in case II when a' = 7'. Thus, quite remarkably, in these cases the 
deviation of | < m > | from the minimal value it can have in the case of neutrino mass spectrum 
with inverted hierarchy (see, e.g., [21]), min(|<m>|) ~ mcos 26*501, is determined by Jcp- 

A similar analysis can be performed and analogous results can be obtained if the small term 
~ e, generating the requisite splitting between the two non-zero eigenvalues of the matrix (|43|) . 
is present in one of the diagonal entries in Eq. (|43|) . 



5 Conclusions 

Assuming that the three flavour neutrino mixing matrix t/pMNS a product of two unitary 
matrices U\ep and 11^ arising from the diagonalization of the charged lepton and neutrino mass 
matrices, ?7pmns = Uj^pUu, and that apart from possible CP-violating phases, Uiy has a bi- 
maximal mixing form, Eq. ((T)), we have performed in the present article a systematic study of 
the possible forms of the matrix Ul^^ which are compatible with the existing data on neutrino 
mixing and oscillations. The three mixing angles in f/iep were treated as free parameters with- 
out assuming any hierarchy relation between them. The case of CP-nonconservation was of 
primary interest and was analyzed in detail. 

We have found that there exist three possible generic structures of ?7iep, which are compatible 
with the existing data on neutrino mixing. One corresponds to a hierarchical "CKM-like" 
matrix. In this case relatively large values of the solar neutrino mixing angle ^soi, and of \Uez\^ = 
|(C/pMNs)e3p, are typically predicted, tan^ ^soi ~ 0.42, |f/e3p ^ 0.02, while the atmospheric 
neutrino mixing angle 6'atm can deviate noticeably from 7r/4, sin^ 26'atm ~ 0.95. The second 
corresponds to one of the mixing angles in U\ep being equal to tt/2, and predicts practically 
maximal atmospheric neutrino mixing sin^ 2^atm — 1- Large atmospheric neutrino mixing, 
sin^ 26'atm ~ 0.95, is naturally predicted by the third possible generic structure of f/iep, which 
corresponds to all three mixing angles in U\ep being large. 

The parametrization of the 3 flavour neutrino mixing matrix as t/pMNS = U^^pU^, is useful 
especially for studying the case of CP-nonconservation. We have found that, in general, in 
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addition to the three Euler-hke angles, f/j^p contains only one CP-violating phase, while f/j^ 
includes four CP-violating phases and can be parametrized as follows: Uu = Pu f^imax Qf, 
where P„ and Q,, are diagonal phase matrices containing two phases each, and f/bimax is the 
bimaximal mixing matrix, Eq. (^. The two phases in Q^, are Majorana-like — the flavour 
neutrino oscillation probabilities do not depend on them. In general, the three angles and the 
Dirac CP-violating phase in the PMNS matrix t/pMNS depend in a complicated manner on the 
three angles and the phase in Uiep and on the two phases in P^. Finally, the two Major ana 
CP- violating phases in f/pMNS are functions of the all the five CP-violating phases in Uiep, Pu 
and Qu, and of the the three angles in Uiep. 

We have derived approximate and rather simple expressions for the observables tan^ ^soi, 
lUesl and sin^ 26'atm, as well as for three rephasing invariants of the ?7pmns matrix, Jcp and S'1,2, 
associated respectively with CP-violation due to the Dirac phase and due to the two Majorana 
phases in t/pMNS- This is done for each of the three possible generic structures of Uiep- In certain 
specific cases we find that simple relations between Jcp and hold (see, e.g., Eq. ^E^). 

Finally, we considered a simple model in which t/pMNS is expressed in terms Ul^^Uy with 
Uy = Pu f/bimax Qv In this model the neutrino mass term is assumed to have the same form 
(when we set Ui^p = 1) as that of the Zee model. One finds that in this case Qv = 1, while 
Py contains two CP-violating phases. Expressions for the earlier indicated neutrino mixing 
observables and for the effective Majorana mass in (/3/3)o,y-decay, |<m>|, are given. We find 
that if the CP-violating phases in P^ satisfy a certain condition, |<m>| depends in a simple 
way on Jcp (Eq. (jSH)). 

It would be interesting to generalize to the matrix Uicp Pu the idea of unitarity triangles, 
which has been successfully used in analyzing CP-violation in the quark sector, particularly 
in B-meson decays. This could provide a helpful graphical representation of CP-violation in 
neutrino oscillations. 
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Figure 1: Scatter plot of the three A parameters for the Sex (left column) and la (right column) 
allowed ranges of values of the neutrino mixing parameters given in Eq. (0)). Conservation of 
CP is assumed. 
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Figure 2: Scatter plot of the A and neutrino oscillation parameters in the case of "small" 
Xij for the la allowed ranges of values of the neutrino mixing parameters given in Eq. (^. 
Conservation of CP is assumed. 
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Figure 3: Scatter plot of the two A parameters and the oscillation parameters in the case of 
A23 = 1 for the la range of the neutrino mixing parameters given in Eq. Conservation of 
CP is assumed. 
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Figure 4: Scatter plot of the three A parameters and the CP-violating phases for the la allowed 
ranges of values of the neutrino mixing parameters given in Eq. (0]). 




Figure 5: Scatter plot of the oscillation parameters for the case of hierarchical Xij, i.e., A13 = 
A12 A23 = Afj- The various possible regions and correlations are shown both in the case of CP- 
violation (bright green areas) and CP-conservation (dark red areas). The la allowed ranges of 
values of the neutrino mixing parameters given in Eq. Q were used. 
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Figure 6: Scatter plot of the correlations between tan^ 6'soi and |f/e3p, and between Jcp and 
|C/e3p; when the Ajj are of the same order and some of the CP-violating phases, ip, and tu, are 
zero. In the last row we display the CP-violating cases corresponding to u; = (bright green 



